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Since the announcement of the proton spin crisis by the European Muon Collaboration 
there has been considerable progress in unravelling the distribution of spin and orbital 
angular momentum within the proton. We review the current status of the problem, 
showing that not only have strong upper limits have been placed on the amount of po- 
larized glue in the proton but that the experimental determination of the spin content 
has become much more precise. It is now clear that the origin of the discrepancy be- 
tween experiment and the naive expectation of the fraction of spin carried by the quarks 
and anti-quarks in the proton lies in the non-perturbative structure of the proton. We 
explain how the features expected in a modern, relativistic and chirally symmetric de- 
scription of nuclcon structure naturally explain the current data. The consequences of 
this explanation for the presence of orbital angular momentum on quarks and gluons is 
reviewed and comparison made with recent results from lattice QCD and experimental 
data. 

1. Introduction 

Amongst the many fundamental questions we face concerning the quark and gluon 
structure of nucleons and nuclei, few problems have proven more exciting than that 
of mapping the distribution of the spin of the proton onto these constituents 
This began in earnest with the announcement by the European Muon Collaboration 
(EMC) — in 1988 that the quarks carried only a very small fraction of the proton 
spin, possibly zero. This soon became known as the "proton spin crisis" . 

Although many effects related to well known features of hadron structure were 
explored none were able to produce such a small result. On the o ther h and, 
it was soon realized that, because of the famous U(l) axial anomaly l ^ l ? ! ^ po- 
lariz ed gluons could also make a contribution to the proton spin structure func- 
tion I 9 | 10 | H | Q 1 Indeed, it was shown that if the gluons in a polarized proton were 
to carry 4 units of spin, the EMC spin crisis would be resolved. No reasonable 
explanation of how the gluons might carry 4 units of spin was ever given. Neverthe- 
less, the mathematical beauty of the proposed resolution was seductive and major 
experimental programs were launched to investigate it. 

Over the past twenty years that experimental effort has produced spectacular 
results, with remarkable increases in precision and new kinds of data which probe 
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for the presence of polarized gluons directly. There have also been some remarkable 
advances in our ability to solve QCD on a lattice of points in space-time and extract 
new information about nucleon structure. Indeed, the progress made has reached 
the stage where it is possible to assert that the gross features of the proton spin 
problem are now understood. This does not mean that considerable experimental 
and theoretical effort will not be required to confirm the resolution. However, as 
we shall explain, well known and understood aspects of hadron structure, in com- 
bination with the totality of modern experimental evidence leave little room for 
doubt. 

We begin by reviewing the experimental situation for the proton spin structure 
function as well as the search for polarized glue - all in the context of the explanation 
proposed in terms of the contribution to gi(x, Q 2 ) through the axial anomaly. We 
then review the possible explanations in terms of more conventional aspects of 
hadron structure which were suggested within a few months of the EMC preprint. 
These explanations are then examined in the light of modern insights into hadron 
structure coming from the study of their properties as a function of quark mass using 
lattice QCD. After summarizing the resolution of the spin problem, we examine the 
consequences for the internal structure of the proton and how this might be tested 
against both experiment and lattice QCD. 

2. Summary of the current experimental situation 

The major advance made by EMC was to extend the measurements of the proton 
spin structure function, g\(x,Q 2 ), to a value of Bjorken-x almost as low as 10~ 2 . 
This enabled a more accurate evaluation of the integral of the spin structure func- 
tion. The latter can be rigorously expressed in terms of well known, perturbative 
QCD corrections and three low energy constants, the axial charges g^ and g^\ 
which are known from neutron and hyperon /3-decay, as well as the quantity S. 
(This quantity, which is often denoted AS or a in the literature, represents the 
renormalizqation group invariant spin content, T,{ji 2 = oo).) Within the parton 
model £ = Au + Ad + As, with Aq — Aq(x)dx, the fraction of the helicity of 
the proton carried by the quarks and anti-quarks of flavor q. (We stress that while 
the word spin can be ambiguous, in the context of the spin problem it rigorously 
refers to the helicity of quarks in a proton, in the infinite momentum frame, with 
positive helicity.) 

In the interest of clarity, we denote the perturbative QCD coefficients, which are 
power series expansions in the strong coupling constant, a s (Q 2 ), as Ci(Q 2 ). Then 
the sum rule for the proton spin structure function is: 



Since g A ' are well known, having an accurate experimental determination of the 
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left-hand side allows one to deduce E. The result of the EMC measurement was 

E = 14 ± 3 ± 10% , (2) 

which is clearly incompatible with 100%, or even 65%, which was the more realistic 
target after accounting for t he r elativistic moti on o f the quarks expected, for exam- 
ple, from bag model studies 021 - see also Ref. 031 f or a similar result in a modern, 
relativistic confining model. Psychologically, it was even more important that this 
result was consistent with zero and this led to considerable excitement concerning 
the Skyrme model, which was initially believed to make just such a prediction. We 
note that it is now known that this is not the case and furthermore, as we shall see, 
the experimental data is no longer compatible with zero. 

In modern terms, Eq. {T]) is especially interesting as one can rigorously extract 
E from a measurement of the axial weak, n eutr al current charge of the proton after 
correcting for heavy quark radiative effects US In this way it becomes a very inter- 
esting sum rule which can really only be violated if there is (5-function contribution 
at x = - as proposed by Bass in connection with the fascinating idea that 
some of the spin of a constituent quark might be associated with topological charge 
stored in its gluon fields. 



2.1. Role of the U(l) axial anomaly 

Although the gluon carries no electric charge, it can create a quark-anti-quark pair 
which then absorbs the incoming photon. If we view the deep-inelastic scatter- 
ing cross section as the imaginary part of the forward photon-nucleon Compton 
amplitude, then this gluonic term corresponds to the imaginary part of a quark- 
anti-quark box diagram. The moments of the structure functions correspond to 
integrating over the momentum of the quark (or anti-quark) which is struck by 
the incoming photon and this effectively shrinks the box diagram to a triangle. 
In the case of spin dependent deep-inelastic scattering, the corresponding integral 
is quadratically divergent and the answer depends on whether or not the regular- 
izati on p rocedure respects gauge invariance or chiral symmetry; it cannot respect 

bothECl. 

If one chooses to work within a scheme that ensures chiral symmetry, the 
renormalized quark spin operators satisfy the usual SU(2) commutation relations 
([Si,Sj] — etjk Sk) and the axial current is conserved. On the other hand, if one 
chooses to ensure the more fundamental symmetry, namely gauge symmetry, the 
renormalized spin operators do not satisfy the usual co mmut ation relations and the 
divergence of the renormalized axial current is not zero HSlSl. In that case, the value 
of E, defined through the matrix element of the axial charge, renormalized in a 
gauge invariant manner, can be written: 

s = s naiV c- Ar/a ; (Q2) AG(g 2 ). (3) 

Z7T 
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Here £ na ive is the spin carried by the quarks in a naive quark model (i.e. one that 
respects chiral symmetry and ignores the effect of the axial anomaly) and the second 
term is the contribution arising, through the axial anomaly, from polarized gluons 
in the proton. 

On the face of it, as a s (Q 2 ) vanishes as Q 2 approaches oo, it appears that the 
correction term on the right of Eq.Q should vanish in the Bjorken limit. However, 
again because of the axial anomaly, the evolution of AG(Q 2 ) is such that the 
product, a s (Q 2 ) AG(Q 2 ), should go to a non-zero constant as Q 2 — > oo. Thus the 
gluon box diagram must be included in the Bjorken limit. 

At a typical scale of Q 2 = 3GeV 2 , the strong coupling, a s (Q 2 ), is around 0.3. 
Thus for 3 active flavors of quark the correction to the naive spin content is roughly 
0.15 AG and if AG were of order 4, the naive expectation for the fraction of the 
spin of the proton carried by its quarks, namely about 2/3 after allowance for the 
relativistic motion of the valence quarks, would be brought into agreement with the 
EMC data. 

It is important to realize a couple of things. First, no-one ever proposed a model 
of nucleon structure in which one could see how the gluons might carry so much 
spin. Secondly, it was understood that if the gluons carried 4 units of spin they 
must also carry about 4 units of orbital angular momentum so that the total spin 
of the nucleon would remain at 1/2. This second caution has tended to be forgotten 
as the limits on the amount of spin carried by gluons have decreased. 

Finally, we note that in the case of an explicitly gauge invariant scheme, such as 
MS, there is no explicit gluon correction to S. Instead, the individual values of Aq 
incorporate the effect of the axial anomaly. Whether one one works in a scheme that 
preserves chiral symmetry or gauge symmetry, provided the effect of the anomaly 
is treated consistently, the same physics is incorporated. It just appears in different 
places. 

2.2. Shape of the gluon contribution to g\ 

Although most of the focus on the gluon contribution to the spin structure function 
of the proton has involved its first moment, it is also of considerable importance 
to know the characteristic shape of the cont ribu tion from the photon-gluon box 
diagram to the proton spin structure function In fact, as emphasized very early 
by Bass and co-workers, the shape of this contribution is distinctive, with 8g\(x) 
going rapidly large and negative^ for x below 10 -2 - if AG is positive, as originally 
required to explain the spin crisis. 

We also note that £ is a flavor singlet quantity, so it will be the same for the 
proton and the neutron. Thus, while the large iso-vector contribution (correspond- 
ing to g A in Eq. {1} ) cancels in the measurement of the deuteron spin structure 
function, the fractional contribution of the photon-gluon box diagram is enhanced. 
This makes the recent, very precise measurements of the deuter on spin structure 
function down to below 10~ 4 , by the COMPASS Collaboration UM^ } especially 
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significant. Those measurements are completely consistent with zero for x from 
5 x 10~ 5 to almost 10 _1 . This alone puts a very severe constraint on the experi- 
mental value of AG - making it impossible for AG to be anywhere near as large 
as 4. 

At Jefferson Lab, the CLAS Collaboration has more than doubled the number 
of polarized spin structure function data in the deep-inelastic region. Because of the 
relatively low range of Q 2 it is cruci al t o allow for higher-twist effects in amalyzing 

loni ioi I 

the data. Nevertheless, Leader et aL'^'have found that including the JLab data 1 ^ 
in their global analysis reduces the error on the determination of xAG by a factor 
of 3 from what is possible without it. Their conclusion is that |AG| < 0.3. 



2.3. Direct searches for AG 

Another characteristic of the gluon box diagram is that it tends to be domin ated 
by quark-anti-quark pairs with transverse momentum equal to or larger than Q 
Thus, one can search directly for the existence of polarized gluons by making semi- 
inclusive measurements that explicitly look for high-p^ hadrons in a deep-inelastic 
event. 

Such di rect m easurements have been extensively explored at Hermes ^2 anc i 

19^1941 

COMPASS l ZJ l z ^ l, with negative results. For two different methods of analysis, Her- 
mes has found AG/G to be of order 0.071 ± 0.034 ± 0.011 and 0.071 ± 0.034 ± 0.010 
in the low- a; region. This led them to the conclusion that AG is small and 
"unlikely to solve the puzzle of the missing nucleon spin". The COMPASS re- 
sult is completely consistent with Hermes but even more precise, with the re- 
sult AG/G = 0.06 ± 0.31 ± 0.06 at x = 0.13 and [i 2 = 3GeV 2 , using data with 
Q 2 > lGeV 2 and looking at high-p^ hadrons. 

At RHIC both PHENIxESESland STAR 12U2S1 have measured asymmetries in 
polarized proton-proton collisions, either looking for asymmetries (i.e. differences 
in the cross sections for helicities anti-parallel and parallel) in either jet or tt° 
production. The preliminary analysis of the latest PHENIX data prefers a value of 
AG between -0.5 and zero. For STAR the preliminary analysis of Run 6 yielded a 
limit for |AG| below 0.3 (at 90% confidence level) and again consistent with zero. 



2.4. Latest status of the spin sum rule 

The tremendous progress in the measurement of both the proton and neutron spin 
structure functions over the twenty years since the original EMC discovery have 
led to a significant increase in the accuracy with which the integral of g\ can 
be determined. We have already mentioned that the range of x over which data 
exists now extends below 10 -4 . In addition, in comparison with the Spin Muon 
Collaboration, the successor to EMC at CERN which was designed to follow up 
the discovery of the spin crisis, the latest data is as much as an order of magnitude 
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more precise. The latest analyses of Hermes ^ and COMPASS ^ yield values: 

£ = 0.330 ± O.Oll(thry) ± 0.025(exp) ± 0.028(evol) Hermes 

S = 0.33 ± 0.03(stat) ± 0.05(syst) COMPASS . (4) 

This represents a very substantial increase in the fraction of the spin of the proton 
carried by its quarks and anti-quarks. 

2.5. Summary 

The result of the last twenty years of intense experimental effort has been a very 
significant improvement in the accuracy with which £ is known. Contrary to the 
original EMC measurement it is no longer possible that the quarks carry none of the 
spin of the proton. Rather it seems that the quarks and anti-quarks in the proton 
carry about a third of its spin and possibly as much as 40%. 

As well as this shift in the target which theory needs to explain, the experimental 
effort has also established very strong upper limits on the spin of the gluons. It 
seems that AG cannot be larger than 0.3 and may even be slightly negative. This 
value of AG is more than an order of magnitude larger than the value (namely 4) 
originally found necessary to reduce the naive theoretical expectation of 65% (after 
allowing for relativistic motion of the valence quarks) to the value observed by 
EMC, through the axial anomaly. Through Eq. ([3]) we see that |AG| < 0.3 implies 
that through the axial anomaly the gluons yield a correction of less than 5% to the 
quark spin content of the proton. Thus the gluon spin is also a factor of 6 too small 
to explain the difference between 65% and the current experimental values of £, 
given in Eq. ((4|). Indeed, if the hint from PHENIX is right and AG is negative, it 
may actually make the problem marginally worse. 

The nature of the challenge presented by the proton spin structure function has 
changed but it is still a fascinating problem of the highest importance. 

3. Modern theoretical explanation 

The discussion in this section has its origins in two papers written in 1988. The 
first, written with Andreas Schreiber SI, dealt with the effect of the pion cloud 
of the nucleon, which is well known to be associated with chiral symmetry. The 
second, written with Fred Myhrer^l, dealt with what in nuclear physics would be 
described as an exchange current correction - in this case arising through the hyper- 
fine one-gluon-exchange (OGE) interaction, which is usually considered primarily 
responsible for the difference in mass between the nucleon and A, the £ and A 
and so on. We consider these two terms in order, noting that our discussion follows 
closely that of Myhrer and Thomas ^H. 

3.1. The pion cloud 

That virtual pion emission and absorption can play a major role in the properties of 
hadrons has been appreciated since its existence was proposed by Yukawa. Within 
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the modern context of QCD as the fundamental theory of the strong interaction 
this remains true, even though we know that the pion is itself composed of quark- 
anti-quark pairs. Indeed, we understand that the unusually small mass of the pion is 
associated with the very low values of the u and d quark masses, with vanishing 
linearly as the light quark masses go to zero. That is, it is a pseudo-Goldstone boson 
associated with the approximate chiral invariance of QCD. 

This special property of the pion leads to many important consequences. For 
example, in the chiral limit (vanishing light quark masses) the charge radius of 
the proton and neutron go to infinity (plus and minus, respectively). Chiral per- 
turbation theory provides a systematic way of exploring the consequences of chiral 
symmetry. In particular, it has established that hadron properties, such as charge 
radii, magnetic moments and so on, are non-analytic functions of the light quark 
masses solely because of the contributions of Goldstone boson loops. Tracing the 
non-analytic behaviour of hadron properties has proven a powerful tool in guid- 
ing the development of models and in testing whether apparent relationships be- 
tween phy sical observables are purely accidental or whether they may be somewhat 
deeper 

An unfortunate consequence of the mathematical elegance of the formalism 
of chiral perturbation theory is that it has been allowed to obscure a number of 
instances where there is a simple, physically meaningful explanation of some inter- 
esting physical phenomena. We recently presented arguments as to why it is both 
meaningful and even satisfying to think of hadron stru cture in terms of a pion cloud 
and we refer interested readers to that discussion 1^1 

In fact, describing a physical nucleon as having a pion cloud which interacts with 
the valence quarks of the quark core (the "bare" nucleon), in a manner dictated 
by the requirements of chira l symm etry, has been very successful in de scribing 
the properties of the nucleon EEl] T ^ e c i ouc jy bag model (CBM) ESSH] renec t s 
this description of the nucleon and in this model the nucleon consists of a bare 
nucleon, \N >, with a probability Z ~ 1 — Pn-h — Pa-k ~ 0.7, in addition to being 
described as a nucleon (iV) and a pion and a A and a pion, with probabilities 
P Nt . ~ 0.20 — 0.25 and Patt ~ 0.05 — 0.10, respectively. The phenom enological 
constraints on these probabilities were discussed, for example, in Refs. I 33 | 36 | 37 [ 
One of the most famous of these constraints is associated with the excess of d over 
u quarks in the proton, predicted on the basis of the CBM ^^Hl. Indeed, to first 
order the integral of d{x) — u{x) is 2/3Pjv7r — Patt/3, which is consistent with the 
experimental data^Slif Pntt and Pa?t he within the ranges just quoted. 

The effect of the pion cloud on the spin carried by quarks was investigated by 
Schreiber and Thomas within a few months of the EMC preprint. They wrote the 
corrections to the spin sum-rules for the proton and neutron explicitly in terms of 
the probabilities set out above 4. For the present purpose it is helpful to rewrite the 
results of Ref. ^ for the proton and neutron. In fact, if we consider explicitly the 
flavor singlet combination, proton plus neutron, the pion cloud modifies the quark 
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contribution to the proton spin in the following manner: 

S - \ Z - ^P N , + ^PaJ £ ■ (5) 

The critical feature of the role of the pion cloud shown in Eq. (JSJ) is that the 
Clebsch- Gordon algebra for coupling the spin of the nucleon and the orbital angular 
momentum of the pion in the Nit Fock state favors a spin down nucleon and a pion 
with +1 unit of orbital angular momentum in the z-direction. This has the effect 
of replacing quark spin by quark and anti-quark orbital angular momentum. Note 
that in the Air Fock component the spin of the baryon tends to point up (and 
the pion angular momentum down), thus enhancing the quark spin. Nevertheless, 
the wave function renormalization factor, Z , dominates, yielding a reduction by a 
factor between 0.7 and 0.8 for the range of probabilities quoted above. 



3.2. The one- gluon- exchange hyperfine interaction 

It is well established that the spin-spin interaction between quarks in a baryon, aris- 
ing from the exchange of a single gluon, explains a major part of the ma ss diffe rence 
between the octet and decuplet baryons - e.g., the N-A mass difference [SES. This 
spin-spin interaction must therefore also play a role when an external probe inter- 
acts with the three-quark baryon state. That is, the probe not only senses a single 
quark current but a two-quark current as well. The latter has an intermediate quark 
propagator connecting the probe and the spin-spin interaction vertices, and is sim- 
ilar to the exchange-current corrections which are well known in nuclear physics. 

In the case of spin dependent properties, the probe couples to the various cur- 
rents in the nucleon. In the first exploration of the two-quark current, illustrated in 
Fig. [TJ carried out by Hogaasen and Myhrer SB the MIT bag model was used and 
the quark propagator was written as a sum over quark eigenmodes. The dominant 
contributions were found to com e from the intermediate p-wave anti-quark states. 

The primary focus of Ref.SB was actually the OGE exchange current corrections 
to the magnetic moments and semi-leptonic decays of the baryon octet, where, for 
example, this correction is vital to understand the unusual strength of the decay 
£~ — > n + e~ +P e . Myhrer and Thomas ^realized the importance of this correction 
to the flavor singlet axial charge and hence to the proton spin, finding that it reduced 
the fraction of the spin of the nucleon carried by quarks, calculated in the naive 
bag model, by 0.15 - i.e., S — * £ — 3G^. The correction term, G, is proportional 
to a s times certain bag model matrix elements SB where a s is determined by the 
"bare" nucleon-A mass difference. 

As in the case of the correction for relativistic motion of the valence quarks 
and for the pion cloud, here too the spin lost by the quarks is compensated by 
orbital angular momentum - here it is orbital angular momentum of the quarks 
and anti-quarks (the latter predominantly u in the p-wave). 
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Fig. 1. The one-gluon-exchange correction to the spin sum rule investigated by Myhrer and 
Thomas 



3.3. Summary 

We have examined three corrections to the fraction of the spin of a nucleon carried 
by its quark and anti-quarks that arise naturally in any realistic treatment of proton 
structure: 



The valence quarks must be treated relativistically. For example, in the 
MIT bag model they satisfy the Dirac equation and the lower component 
of the Dirac spinor for a spin-up valence quark in an s-state will have 
predominantly L z = +1 and spin down. This reduces the fraction of the 
spin carried by valence quarks from 100% to around 65%. 
The Clebsch-Gordon algebra for the dominant Fock component of the nu- 
cleon wave function associated with its pion cloud is identical. That is, the 
dominant term has the quark core (or "bare nucleon" , with spin down) 
while the pion has predominantly L z = +1. Within the context of chiral 
quarks models, such as the CBM, this means that whatever the spin con- 
tent within the quark model used, dressing it with the pion cloud reduces 
that value by a factor of order (0.7,0.8). 

The exchange current correction associated with OGE in the proton also 
reduces the amount of spin carried by the quarks by about 15%, where this 
number has so far only been calculated within the MIT bag model. 
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Clearly, if we combine these three corrections, that is we assume a relativistic 
description of the structure of the nucleon which incorporates OGE and chiral 
symmetry, the simple accounting of the fraction of the nucleon spin carried by its 
quarks and anti-quarks is: 

S = (0.7,0.8) x (0.65 - 0.15) 

-> £ e (0.35,0.40). (6) 

This result is in very satisfactory agreement with the current experimental values 
given earlier in Eq. 0$ and the problem appears to be solved - at least at the level 
of a few percent. 

There were two reasons why this conclusion could not be drawn in 1988. First the 
experimental target was 14% and possibly zero and even combined these corrections 
could not explain that. Secondly, at the time, the various studies of chiral quark 
models suggested that a large fraction of the mass splitting between the N and A 
arose from pion cloud corrections. If that were so, it would be double counting to 
include both the pion cloud correction and the OGE correction with a strength 
determined by the observed N-A mass splitting. 

Experimental progress over the past 20 years has eliminated the first problem. 
For the second, we owe a resolution to the sophisticated studies of hadron properties 
as a function of quark mass that have been stimulated by lattice QCD over the past 
decade. As a result of discoveries by Leinweber, Young and Thomas concerning 
the link between quenched and full lattice QCD, we now know that only 40 ± 20 
MeV of the observed 300 MeV N-A mass splitting comes from the pion cloud SH. 
Consequently, there is little or no double counting as 80-90% of the N-A mass 
difference would then come from OGE and the effective strong coupling constant, 
a s , used to calculate the OGE correction to the quark spin is appropriate. 

In the next section we review the relevant information concerning the analysis 
of lattice QCD data as a function of quark mass, which led to the discovery just 
mentioned. 

4. Origin of the N-A mass difference 

One of the unexpected but very positive consequences of our lack of supercomputing 
power is the fact that it has not been possible to compute physical hadron properties 
in lattice QCD. In fact, with computation time scaling as a large inverse power of the 
pion mass, calculations have covered the pion mass range from 0.3 to 1.0 GeV (or 
higher). Far from being a disappointment, this has given us a wealth of unexpected 
insight into how QCD behaves as the light quark masses are varied ™l. In terms of 
the insight this has given us into hadron structure it is both truly invaluable and 
thus far under-utilized. 

The most striking feature of the lattice data is that in the region > 0.4 GeV, 
in fact for almost all of the simulations made so far, all baryon properties show a 
smooth dependence on quark mass, totally consistent with that expected within a 
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constituent quark model. The rapid, non-linear dependence on m n required by the 
LNA and NLNA behavior of xPT are notably absent from the data! 

The conventional view of xPT has no explanation for this simple, universal 
observation. Worse, in seeking to apply x?T to extrapolate the data back to the 
physical pion mass, it has been necessary to rely on ad hoc cancellations between 
the high order terms in the usual power series expansion (supplemented by the 
required non-analytic behavior). In fact, there is strong evidence that such s eries 
expansions have been applied well beyond their region of convergence I 44 | 45 | an( j 
that, as a result, the corresponding extrapolations are largely unreliable. 

On the other hand, the picture of the pion cloud that we have discussed briefly 
yields a very natural explanation of the universal, constituent quark model be- 
haviour of hadron properties found in the lattice simulations for ra, > 0.4 GeV. 
The natural high momentum cut-off on the momentum of the emitted pion, which 
is associated with the finite size (typically R ~ 1 fm) of the bare baryon (i.e., the 
bag in the CBM), strongly suppresses pion loop contributions as increases. The 
natural mass scale which sets the boundary between rapid chiral variation and con- 
stituent quark type behavior is 1/R ~ 0.2 to 0.4 GeV. Indeed, when in the early 
investigation of the quark mass dependence of nucleon properties the CBM was 
compared directly with lattice data, the agreement was remarkably good S3. (Sim- 
ilar results have been obtained recently within the chiral quark soliton model SSI ) 
The r esults w ere equally as impressive for the N and A masses and magnetic mo- 
ments ^®|49 l the proton charge radius and the moments of its parton distribution 
functions The key features necessary to reproduce the behaviour found at large 
quark mass in lattice QCD and to reproduce the experimentally measured data at 
the physical mass seem to be that: 

• The treatment of the pion cloud (chiral) corrections ensures the correct 
LNA ( and NLNA, although in practice this seems less important in many 
applications) behaviour of QCD 

• The pion cloud contribution is suppressed for beyond 0.4 GeV, and 

• the underlying quark model exhibits constituent quark like behaviour for 
the corresponding range of current quark masses. 

While the CBM satisfies all of these properties, in analyzing lattice QCD data 
one does not want to rely on any particular quark model and this led to early inves- 
tigations of interpolating formulas which built in th e correct asymptotic behavior, 
both in the chiral regime and the larg e mass limit I 51 | 52 | 53 | However, within the 
framework of effective field theory it is more appropriate to suppress the pion 
cloud as m,r goes up by using a simple, finite range regulator (FRR, also refered to 
as a long distance regulator by Donoghue et al. ^j) in the evaluation of the pion 
loops. In this way one can ensure the correct LNA and NLNA behaviour, as well as 
the right large mass behavior, at the cost of one additional parameter, the cut-off 
mass A. If the data are good enough one can use this as a fitting parameter but in 
general it is sufficient to choose a value consistent with the physical arguments 
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presented above. The sensitivity of the extrapolation to the choice of the functional 
form of the FRR is then an additional source of systematic error in the final quo ted 
result. In the case of the nucleon mass the corresponding systematic error was^^ 
of the order of a mere 0.1%. 

One of the most remarkable results of this physical understanding of the role 
of the pion cloud and, in particular, its suppression at large pion mass has been 
the unexpected discovery of a connec tion between lattice simulations based upon 
quenched QCD (QQCD) and full QCD^. In a study of the quark mass dependence 
of the N and A masses SH, it was discovered that if the self-energies appropriate 
to either QQCD or full QCD were regulated using the same dipole form for the 
FRR (the dipole being the most natural physical choice given that the axial form 
factor of the nucleon has a dipole form) with mass parameter A = 0.8 GeV (the 
preferred value), then the residual expansions for the nucleon mass in QQCD and 
QCD (and also for the A in QQCD and QCD) were the same within the errors 
of the fit! This is a remarkable result which, a posteriori, gives enormous support 
to the physical picture of the baryons consisting of confined valence quarks sur- 
rounded by a perturbative pion cloud. The baryon core is basically determined by 
the confinement mechanism and provided the choice of lattice scale reproduces the 
physically kno wn confining force (either through the string tension or the Sommer 
parameter derived from the heavy quark potential) it makes little difference 
whether one uses QQCD or full QCD to describe that core. What does matter is 
the change in the chiral coefficients as one goes from QQCD to full QCD. 

Perhaps the most significant application of this discovery has been the appli- 
cation to the calculation of the octet magnetic moments and charge radii based 
on accurate QQCD simulations th at extend to rather low quark mass. Using the 
constraints of charge symmetry I^Sl this has led to some extremely accu rate calcu- 
lations of the strange quark contributions to the magnetic moment ^01 and charge 
radius of the proto n. Inde ed, those calculations are in excellent agreement with 
the current world data ^^2l but, in a unique example in modern strong interaction 
physics, they are an order of magnitude more accurate. 

In the present context, the key result of the analysis of QQCD and full QCD 
data for the N and A masses is that the contribution of the pion cloud to their 
mass difference is really under control. The result is that only 40 ± 20 MeV of the 
observed mass difference can be attributed to pions. The rest must be associated 
with other mechanisms such as the traditional OGE hypefine interaction. 



5. Quark orbital angular momentum 

The unifying feature of each of the mechanisms which we have discussed, leading to 
the reduction of the spin sum rule, is that quark spin is always replaced by quark 
or anti-quark orbital angular momentum I 31 | 64 |_ -j,^ success i ve effect within the 
proposal of Myhrer and Thomas, as one goes from a non-relativistic quark picture, 
to a relativistic quark model and then includes OGE and pion loops, is illustrated in 
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Table [T] We see clearly that the spin of the proton resides predominantly as orbital 
angular momentum of the u (and u) quarks. In contrast, the d (and d) quarks 
carry essentially no orbital angular momentum. The total angular momentum is 
shared between the u (and u) quarks, J", and the d (and d) quarks, J d , in the 
ratio J" : J d — 0.74 : —0.24. (Note that there are no strange quarks in the Myhrer- 
Thomas calculation, so £ in Table 1 is Au + Ad. Combining this with g\ = Au — 
Ad = 1.27 yields these values. A more sophisticated treatment, including the KN 
Fock compo nen t of the proton wavefunction^Sl, would lead to a very small non-zero 
value of As^.) 



Table 1. Distribution of the fraction of the spin of the nucleon as the spin and orbital angular 
momentum of its constituent quarks (and anti-quarks summed) at the model (low energy) scale. 
Successive lines down the table show the result of adding a new effect to all the preceding effects. 



L u L d S 

Non-relativistic 1.00 

Relativistic 0.46 -0.11 0.65 

OGE 0.67 -0.16 0.49 

Pion cloud 0.64 -0.03 0.39 



5.1. Comparison with lattice QCD 

Over the past few years there have been extensive studies of the low moments of the 
so-called Generalized Parton Distributions, which are directly related to the total 
angular momentum carried by the quarks within the nucleon I 67 | 68 [ Those stud- 
ies suggest that L d tends to be positive, while L u is negative. One should observe 
that these calculations were performed at fairly large quark mass and omit discon- 
nected terms, which may carry significant orbital angular momentum and are 
certainly needed to account for the U(l) axial anomaly. Nevertheless, the apparent 
discrepancy is of concern. 

In order to resolve this problem, as pointed out in Ref. ^1, it is important to 
recall that neither the total, nor the orbital angular momentum is renormalization 
group invariant (RGI) ^1 The lattice QCD values are evaluated at a scale, set by 
the lattice spacing, around 4 GeV 2 . On the other hand, we have not identified the 
scale corresponding to the values derived in our chiral quark model. Indeed, there 
is no unambiguous way to do so unless the model can be derived rigorously from 
non-perturbative QCD. 

This problem has been considered for more than 30 years driven initially by 
the fact that in a typical, valence dominated quark model, the fraction of momentum 
carried by the valence quarks is near 100%, whereas at 4 GeV 2 the experimentally 
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Fig. 2. Evolution of the total angular momentum and the orbital angular momentum of the up 
and down quarks in the proton, as given at the low starting scale in Table 1 (from top to bottom 
(at 4 GeV 2- !: J " (solid), L d (smallest dashes), L u (largest dashes) and J d (middle length dashes)) 
-fromRef.ES. 



measured fraction is nearer 35%. Given that QCD evolution implies that the mo- 
mentum carried by valence quarks is a monotonically decreasing function of the 
scale, the only place to match a quark model to QCD is at a low scale, Qq. Early 
studies within the bag model found this scale to be considerably less than 1 GeVL^I. 

Over the last decade, this idea has been used with remarkable success to describe 
the data from HERA, over an enormous range of x and Q 2 , starting from a valence 
dominated set of input parton distributions at a scale of order 0.4 GeV^l A similar 
scale is needed to match parton distributions calculated in various modern quark 
models to experimental data^^l. Indeed, one may view the choice of starting scale as 
part of the definition of the model. We note that the comparison between theory and 
experiment after QCD evolution is not very sensitive to the order of perturbation 
theory at which one wor ks. However, what does change is the unphysical starting 
scale. Following Ref. ^1, we show results only at leading order, which also avoids 
questions of scheme dependence. 

The QCD evolution equations for an gula r momentum in the flavor singlet case 
were studied by Ji, Tang and Hoodbhoy^l The scheme used corresponds to the 
choice of a renor mal ization scheme which preserves chiral symmetry, rather than 
gauge symmetry so that £ is scale invariant. In Fig. [2] we show the effect of 
evolution from a model scale of 0.4 GeV to the scale of the lattice QCD data on 
the results reported in Table 1. 

In contrast with the behaviour of J u ' d in Fig. [21 which is is unremarkable, the 



April 10, 2009 19:24 WSPC/INSTRUCTION FILE paper 



Spin and orbital angular momentum in the proton 15 

corresponding behaviour of L u,d is spectacular. L u is large and positive and L d very 
small and negative at the model scale but they very rapidly cross and settle down 
inverted above 1 GeV 2 ! The reason for this behaviour is easily understood, because 
asymptotically L u and L d tend to 0.06- Aw/2 and 0.06- Ad/2, or -0.36 and +0.28, 
respectively. This is a model independent result and it is simply a matter of how 
fast QCD evolution takes one from the familiar physics at the model scale to the 
asymptotic limit. 

Clearly, in spite of the criticisms one could make of the state-of-the-art lattice 
simulations, especially the absence of disconnected terms and the relatively large 
quark masses involved, the results just reported are consistent with the latest lattice 
results of Hagler et al. For example, they report J u+d in the range 0.25 to 0.29 
at the physical pion mass (their Fig. 47) in comparison with 0.30 in the calculation 
reported above. They also report L u+d ~ 0.06 in comparison with 0.11 in this work. 
Finally, the qualitative feature that L d is positive and bigger than L u is, as we have 
explained, clearly reproduced by the Myhrer-Thomas work after QCD evolution. 

5.2. Comparison with experimental data 

The extraction of information about the quark angular momentum from experi- 
mental data on GPDs is still in its very early stage of development. One needs to 
rely on a model to analyze the experimental data, which are still at sufficiently low 
Q 2 that one cannot be sure that the handbag mechanism really do minates . Never- 
theless, the combination of DVCS data on the proton from Hermes l ^ | 75 | an( j -j- ne 
neutron from JLab^^ (both at a scale Q 2 ~ 2GeV 2 ), provides two constraints on 
J u and J d , within the model of Goeke et al. ag snown m Fig. [3j Also shown 

there is the prediction based on the model of Myhrer and Thomas, as calcvulated 
in Ref. ^1 Note that the error bands are the purely experimental (predominantly 
statistical) errors and there is, as yet, no information on the possible systematic 
variation corresponding to a change of model. The exploration of the model depen- 
dence is clearly a high priority for future work. Nevertheless, within the present 
uncertainties, most notably the relatively low Q 2 of the data and the unknown 
model the relatively low Q 2 of the data and the unknown model dependence of the 
extraction of J"W, there is a remarkable degree of agreement. 

5.3. Observation on the gluonic correction to spin sum rule 

Under QCD evolution from the model scale, AG(Q 2 ) grows from zero to a quantity 
that increases linearly with InQ 2 in the asymptotic region. With the starting values 
given in Table 1, AG grows to around 0.5 at Q 2 — 4 GeV 2 . Through Eq. {3]), 
this means that the chiral value of S, calculated by Myhrer and Thomas, S G 
(0.35,0.40), should receive a finite correction before it can be compared with the 
experimentally determined value, £(Q 2 = oo) = Sinv Making this correction using 
leading order QCD evolution leads to a theoretical value for £j nv £ (0.25,0.29). 
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Fig. 3. Comparison between the constraints on the total angular m oment um carrie d by u and d 
quarks in the proton, derived from experiments on DVCS at Hermes EHE] and JLab EST a t a s cale 
of order 2 GeV 2 , and the model of Myhrer and Thomas (the small dark rectangle)- from Ref. EH 

This also is completely consistent with the current experimentally allowed range, 
given in Eq. |4]). 

6. Summary and outlook 

Two decades of experimental investigation have given us a wealth of important 
new information concerning the spin structure of the proton. We now know that 
the spin crisis is nowhere near as severe as once thought but still only about a third 
of the proton spin is carried by its quarks. Polarized gluons, which in principle could 
contribute through the axial anomaly, in practice seem to play no significant role. 
It seems likely that less than 5% of the missing spin can come from polarized glue 
and the sign may be such that it makes the problem slightly worse. 
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Instead it appears that important aspects of the non-perturbative structure of 
the nucleon do resolve the crisis. Indeed, three pieces of physics present in any re- 
alistic model are required. These consist of the relativistic motion of the valence 
quarks, the one-gluon-exchange interaction needed to describe the hyperfme split- 
ting of hadron masses (especially the N-A mass difference) and finally the inclusion 
of the pion cloud required by chiral symmetry. These three terms reduce the fraction 
of the spin of the proton carried by its quarks to between 35 and 40%, in excellent 
agreement with the latest experimental data. If we allow for the small correction 
arising through the axial anomaly from the small amount of radiatively generated 
polarized glue, as explained in Sect. 15.31 the theoretical prediction for the invariant 
spin fraction is in the range (0.25,0.29). The theoretical consistency of this picture 
owes a great deal to recent studies of the dependence on quark mass of hadron 
properties calculated in lattice QCD. 

After QCD evolution from the relatively low scale that typically characterizes 
quark models to 4 GeV 2 , there is a surprising qualitative change in the distribution 
of orbital angular momentum, with the up quarks tending to have a small or even 
negative orbital angular momentum and the down quarks a positive value of order 

0. 1 or larger. This result appears to be consistent with those emerging from recent 
lattice QCD simulations. It will be very important to pursue those calculations to 
more realistic quark masses. 

For the future, it will be critical to test that the missing components of the 
proton spin do indeed reside as orbital angular momentum on the quarks and anti- 
quarks, as implied by this theoretical explanation. In this respect the program of 
measurements of Generalized Parton Dis tribution s, especially following the 12 GeV 
Upgrade at Jefferson Lab, will be vital I 79 | 80 | 81 [ we nave discussed, the early 
results from Hermes and JLab at 6 GeV, albeit within a particular model and at 
an uncomfortably low value of Q 2 , do seem to be consistent with the explanation of 
the spin problem proposed by Myhrer and Thomas. This certainly provides great 
encouragement for further work. 
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